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Abstract: Here we study zero temperature quantum phase transition driven 
by the transverse field for random ± J Ising model on chain and square lattice. 
We present some analytical results for one dimension and some numerical results 
for very small square lattice under periodic boundary condition. The numerical 
results are obtained employing exact diagonalization technique following Lanczos 
method. 



Keywords: Quantum phase transition, Spin glasses, Transverse Ising model, 
Diagonalization techniques. 

PACS Nos.: 05.70.Fh, 42.50.L, 75.10.N 

1. Introduction 

The interest in the study of transverse Ising spin glass models was revived in early 
1990 by the discovery of zero-temperature transition in dipolar Ising transverse 
field magnet LiH02.Yx_j.F4 [1]. Proton glasses such as mixture of ferroelectric and 
anti-ferroelectric materials like Rb(i_ x )(NH4)_(H2P)4 [2] also provided earlier use- 
ful realizations of such quantum spin glasses. 

These developments initiated extensive theoretical studies in quantum spin glass 
models. Ising model in transverse field has already been studied extensively in 
this context through analytical approaches using approximate renormalization 
techniques and real space renormalization group method, as well as using nu- 
merical methods like quantum Monte Carlo and exact diagonalization techniques 
[3]. Fairly extensive studies on the quantum spin glass phases have been made 
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so far for transverse field Edward-Anderson model and Sherrington-Kirkpatrick 
model with random exchange distributions [3, 4] and to some extent on quan- 
tum Heisenberg spin glass models [4]. In view of the rigorous developments in 
the study of two dimensional nearest neighbour Ising model with random ±J 
exchange interactions and the precise knowledge of location of the Nishimori line 
[5] in such classical spin glass model (driven by temperature), we consider here 
the quantum phase transition (at zero temperature) in the same ± J Ising model. 

We have shown analytically that introduction of random — J impurities cannot 
affect the zero temperature phase transition in one dimensional system as they 
can be transformed away. We have also compared and verified the result nu- 
merically for the small system size considered. For two dimensional systems, we 
present some priliminary results obtained for a square lattice using exact diag- 
onalization results for very small system sizes following Lanczos technique [6]. 
Only the behaviours of configurationally averaged energy gap A = (Ei — E ) be- 
tween the first excited state and the ground state and the second order response 
function x — (d 2 E /dT 2 ), equivalent to specific heat, have been studied here. 
The variations of A and \ with respect to transverse field T have been obtained, 
and the phase boudary has been estimated from these results. 

We work with a transverse Ising system, using only nearest neighbour interac- 
tions, whose Hamiltonian is given by 

JV 

w = -EWs/- r Es?> (i) 

where the transverse field T is uniform through out the system and the nearest 
neighbour exchange constant J^'s are chosen randomly from the binary distribu- 
tion 

P(J tJ ) = P 5(J i3 + J) + (1 - p)S(J ij - J). (2) 

Here J is taken positive and p is thus the concentration of anti-ferromagnetic — J 
bonds in the system. 

2. Results in One Dimensional System 

Here first we show analytically that in a one dimensional transverse Ising Hamil- 
tonian with uniform J and T, if we replace some J bonds by — J bonds randomly, 
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then the resulting Hamiltonian can be gauge transformed back to one with uni- 
form J, and hence the critical field remains unchanged with randomness concen- 
tration. Simillar result for one dimensional system with distributed J had been 
obtained earlier [7]. 

Let us take the one dimensional random bond tranverse Ising Hamiltonian 

w = -EW^i-E r <s?> (3) 

i 

where the transverse field T is uniform throughout the system, and Jj's are ran- 
domly chosen from the same distribution as given in (2). Since the Jj's have same 
magnitude J all through, and their randomness is only in their sign, we may write 
Jj = Jsgn(Jj) , and thus Hamiltonian (3) takes the form 

n = - J £ s g n( Jjsts^ - r £ Sf. (4) 

Now let us define a new set of spin variables as below 

St = St ff sgn( J k ) 
k=i 

QX QX 

o i — o i 

Sf = sfllsgn(J k ). 

k=i 

It is easy to see that S 's satisfy the same commutation and anti-commutation 
relations as those of S"s and hence will exhibit exactly the same dynamical be- 
haviour. Now, 



QZ QZ QZ QZ 



i-1 

n[sgn(J fe )] 2 sgn(J), 
k=i 



or, 



since [sgn(Jfc)] 2 = 1. Thus in terms of new spin variables, Hamiltonian (4) be- 
comes 

w = -j£^i- r Esr- (5) 
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The above Hamiltonian describes the same random system in terms of new vari- 
ables, and yet, as one can see, it has in itself no randomness at all. One can 
use Jordon-Wigner transformation in terms of S"s and see that here also quan- 
tum phase transition occurs only at T > T c (= J) as it occurs in a non-random 
Hamiltonian in S"s. In Fig. 1, we present some data computed for a chain of size 
iV = 9, which shows that the gap A vanishes at T c pa 1 (the field being scaled by 
J). These data for A = E\ — E is obtained from the computed average value of 
Eq and E\, each one averaged over about 10 configurations for p ^ 0. For infinite 
system, A is a linear function of V for V > T c . In our case, linearity is observed 
at high values of T, and T c is determined by backward linear extrapolation from 
the linear region. 
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Figure 1: A numerical estimate of the energy gap A for a chain with iV= 9. 
Phase boundary is obtained from the location of A(r) =0, and is shown in 
the inset. 

From the numerical data in Fig. 1, we see that there is a slight variation of T c 
with p (r c varies between 0.9 and 1.0). This variation can be attributed to the 
very small size of the system. However, it may be noted that with even number 
of —J bonds in the chain, with periodic boundary condition, there is no prob- 
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lem of incommensuration and E (T) or -E'i(r) become strictly identical for such 
values of p. Similarly, in every case of odd number of —J bonds in the chain, 
incommensuration problem always occurs for one spin only, rendering identical 
values (but different from the even — J case) for E (T) and -E'i(r) in all such cases. 



3. Results for Square Lattice 

We consider now the same system (represented by Hamiltonian (1)) on a square 
lattice of size 3x3 with periodic boundary condition. We again calculate the 
ground state and the first excited state energy E and Ei respectievely as func- 
tions of the transverse field T, for different values of p. Each value of E and Ei 
is averaged over at least 10 configurations for each p ^ 0. Apart from A, we also 
calculate x — d 2 E /dT 2 and their variations with T as shown in Figs. 2 and 3 
respectively. 
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Figure 2: A numerical estimate of the configurational avrage of the energy 
gap A for a square lattice of size 3x3. Phase boundary obtained from A(r) = 
is outlined in the inset. 



Our results here are severely constrained by the system size. The value of pure 
ferromagnetic critical field T c (p = 0) is found here to be around 2.2, while the 
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Figure 3: Here variation of x — d 2 Eo/dT 2 with T is shown. The transition 
point occurs at T = T c where x diverges; for finite system one gets only a 
peak in x at T = T c (p). We have outlined the corresponding phase boundary 
in the inset. 



series results [3] or cluster algorithms [8] give the value to be around 3.0. This dis- 
crepancy is attributed to the smallness of our system size (N = 3 2 ). However the 
qualitative behaviour of the order- disorder phase boundary (between ferro/spin 
glass and para) seems to be reasonable: T c (p) decreases with p initially, and then 
increases again as p approaches unity (pure anti-ferromagnet). The use of peri- 
odic boundary condition here (to avoid some numerical errors) also restricts the 
domain features and thereby affects our results. The absence of the knowledge 
of the ground-state wave function (and the correlation functions) in this method 
also forbids us to analyse the structure of the ordered phases. 
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